
function [tau_SB1Rot , xb2c_SB1Rot] = Bcr4bp_EMRot2SB1RotMtx3D(tau_EMRot , xb1c_EMRot , aux)
%
% BCR4BP坐标系转换，EMRot -> SB1Rot，3D，矩阵计算
%
% s：太阳
% e：地球
% m：月球
% c：卫星
% b1：地-月质心
% b2：太阳-b1质心
% 
% 参考：
% 2024_Scheuerle, page160
%
% 时间：2024年06月03日
% 作者：张晨
% 单位：CSU
% 邮箱：chenzhang@csu.ac.cn
%%%%%%%%%%%%%%%%%%%%%%%

N = length(tau_EMRot);

tau_SB1Rot = zeros(size(tau_EMRot));
xb2c_SB1Rot = zeros(size(xb1c_EMRot));

for iLoop = 1 : N

    [tau_SB1Rot(iLoop) , xb2c_SB1Rot(iLoop , :)] = Bcr4bp_EMRot2SB1Rot3D(tau_EMRot(iLoop) , xb1c_EMRot(iLoop , :) , aux);

end

end
